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1 Introduction and Motivation 



Constructing even the lowest state that is the 'vacuum' of a Quantum Field Theory has been so 
notoriously difficult that the conventional Hamiltonian approach was given up altogether, long ago 
in the Fifties. It was overlooked that other forms, particularly Dirac's 'front form' of Hamiltonian 
dynamics might have less severe problems. In fact, concretely realizing the front form with 
periodic boundary conditions one might combine the aspects of a simple vacuum [Q] with a careful 
treatment of the infrared degrees of freedom. The method referred to is Discretized Light-Cone 
Quantization (DLCQ) Q. A review can be found in |^. Potentially, the method is able to reconcile 
the simplistic but otherwise successful constituent quark picture of Feynman's parton model [|[ 
with the low-energy regime of Quantum Chromodynamics (QCD) as the fundamental theory of 
hadrons. But the conversion of this non-perturbative method into a reliable tool for hadronic 
physics is bestowed with many difficulties Q, one of them being the so-called 'zero mode problem'. 

Recently, (2-1-1) dimensional pure SU(2) gauge theory has been examined suppressing transverse 
gluon momentum excitations The dimensionally-reduced theory turns out to be SU(2) gauge 
theory in (H-l) dimensions coupled to adjoint scalar matter. A topological (gauge) mode appears 
coupled to true dynamical Fock modes of the transverse gluon fields. Moreover, a constrained zero 
mode appears which is defined by a linear but still very complicated equation, whose inversion is 
far from being trivial [^, structures that were foreseen by Franke et al. |^] in (3-1-1) dimensions. 
Such types of models, first discussed in DLCQ by [|l^, 0, |l^, 13| but without zero modes and 
assuming only color singlet string states, enable insight into how to overcome the obstacles in the 
full theory. 

Next to interesting vacuum structures ||8|, the dimensionally reduced model of QCD has a 
fascinating excitation spectrum whose structure is difficult to guess, see 0. It is this latter point 
to which we address ourselves in the present work. We do that at the expense of rigor, keeping 
only the Fock space structure of the model to be presented below in detail. Despite the heavy 
truncations the present work has rather interesting aspects, which can serve as paradigms for what 
one expects in full QCD. 
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2 Formulation of the Model 



The model ||7| considered in the sequel is a (1+1) dimensional non-abelian gauge theory covariantly 
coupled to scalar adjoint matter whose Lagrangian density is given by 

£ = TY(-iF"^F„^ + D°$D„$-^g^>$). (1) 

The color electro-magnetic fields are the usual ones in matrix notation, i.e. 

F"^ = a"A^-a^A" + ig[A",A'3] =a"A^-D^A". (2) 

The covariant derivative is the carrier of 'longitudinal gauge invariance' which is not aff'ected 
by adding a mass term. Fixing the gauge 

Tr(r±A+)=0, 5_A+ = 0, and TV/r^A") =0, (3) 





completes the model. With the purpose of setting the notation we recall here in short the most 
important steps of Ref.]^]. 

For SU(2), the above r° are 2x2 matrices related to the Pauli spin matrices ai. They define a 
color helicity basis in terms of which all fields are expanded, i.e. 

= ^0-2 , = -^^((^x ± icTy) , thus $ = T^ips + r+v?+ + T~ip_ . (4) 

We shall be careful to have the Lorentz and color indices consistently raised and lowered, re- 
spectively, when we write down indiviual field components. Since the matrices are traceless and 
hermitean, one deals with nine real-valued operator functions, three for each field. Gauge fixing, i.e. 
A+ = VT^, reduces the problem to six real- valued operator functions and one quantum mechanical 
operator v = v{x^), which one keeps track of in the combination z = vgL/ir. Four out of these seven 
have a non- vanishing conjugate momentum, particularly the momenta canonically conjugate to v, 
(^3, if-, and (^+. Only these are therefore 'independent fields' subject to be quantized. The vector 
potential v (thus z) is quantized like a quantum mechanical variable |jl^ that is like [z,Pz] = i, with 
the momentum pz = li^dj^vj g conjugate to z. With 994. = (^L and imposing periodic boundary 
conditions (DLCQ), the remaining two fields are quantized canonically. It is justified Q to restrict 

to the fundamental modular domain (0 < 2 < 1), in which they are represented by 

1 00 

ip3{x-) = ^+ Y^i^nWn e'^'^f + at wn e+^^f-") , (5) 
V47r V47r;fr"i^ ^ 



and ^^{x-) = f^fb^Ur^e-'^T--+dlvme+'^T--). (6) 
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The single particle operators an, bm and dm are the carriers of the operator structure and obey 
conventional commutator relations like [an,a|^/] = (5n,n' and [bm,bl^'] = [dm-,dl^,] = 5m,m.' with all 
others vanishing except those involving the zero mode oq- The coefficients 

Wn = l/y/n , Um = l/ Vm + C , and Vm = 1/ \/ m - C, , with C = z-^, (7) 

are real and depend on the vector potential v through It is justified to think of the 'a'-, 'b'-, and 
'd'-particles as of photons, electrons, and positrons, except that they are all bosons. Below, use of 
this will be made both in the notation and the diagrammatical representation. 

Neither ag, the zero mode of (/?3, nor the remaining three fields ^3 , AZ and A'^ have a conjugate 
momentum. They cannot be quantized canonically but determine themselves by certain equations 
of motion. The A~ are determined by the three Gauss' equations. In terms of the current 



J^ = -[$,D'3$] (8) 

i 

and its density components J+ = 2Tr (r'^J+) they read as 

- = gj+ , + igvfA- = gj+ , and - (a_ - igvfA' = gj+ . (9) 

The first of them can be solved only if the zero mode Q of the r.h.s 

{4\ ^ ^ dx-J+ix-) , thus Qs = 2L = J2 {dlndm - blbm) , (10) 

m=i 

vanishes. This cannot be satisfied as an operator, but must be used to select out physical states. To 
find them is easy: they must have the same total number of 'b' and 'd' particles. When discussing 
below the spectra we will see, that this condition is different from the naive picture of color singlets. 
Finally, one must express ao in terms of the independent fields. The defining constraint equation 
is linear in oq, for details see Ref. Q. In the solution, oq should be an explicit functional of the 
Fock-space operators, and a function of C and fi, i.e. 



ao = ao 



a: 



n 1 ; '^m i ' ^rn ; dm j C ) • ( ^ ^ ) 



Thus far, only approximate solutions have been constructed 
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3 On the Hamiltonian and its Diagonalization 

Upon evaluation the total (light-cone) momentum becomes a diagonal and thus simple 
Fock-space operator 

P+ = - ^n alan + j- ^ [{m + () bl,bm + (m - C) dLdm ■ (12) 

n=l m=l 

The total (light-cone) energy P^ , the 'Hamiltonian' of the theory, is a off-diagonal Fock-space 
operator M. It is so complicated that we only outline its construction and refer for details to M. 



The defining equation is with g = g/^/lOn 

L (f 1 

P- = -Ag^-^ + - jdx-[^il{(t>3(l)3 + (l)+(t>- + (t>-^+) + d-A^d-A:^ 

-L 

+L 

+ 2 ((5_ + igv)Al (5_ - igv)AZ + {d- - igv)AZ (9- +igv)A^y (13) 



-L 



To get it as a Fock-space operator one expresses the charge densities of Eq. (^) in terms of the (pa 
and substitutes Eqs. (^) and (|6|). Inverting Eqs.(^) yields A~ in terms of the Fock-space operators 
and of ao, thus Q 

; ajj, 6j„, dj^, a^) ^mj '^mj C) A* • 0-^) 



P- 



Finally, one has to invert the constraint equation for ao and to substitute Eq.(|ll|) everywhere. As 
a results one gets the Hamiltonian in terms of 'raw-ordered' Fock-space operator products. It is 
reasonable to rewrite these 'time ordered products' as normal ordered products plus the sum of 
all pairwise contractions. The totally contracted terms are identical with the Fock-space vacuum 
expectation value which we conveniently define as W^(C) = (0| -P~ |0)- Contrary to conventional 
DLCQ these may not be discarded since as functions of ( they are part of the operator structure. 
Terms with one creation and one destruction operator are usually referred to as the 'contraction 
terms' or the 'self-induced inertias' Pq , those with two creation and two destruction operators 
as 'seagulls' Pg , those with one creation and three destruction operators as 'forks' Pp. Explicit 
formulas can be found in Appendices ^ and ^ respectively. They have been derived for oq = 0, 
but one should emphasize that a non-vanishing oq contributes to the 'gauge potential energy' T^(C) 
as well as to the contractions, seagulls and forks. In general, the Hamiltonian becomes thus 

TTdC 



P~ = + ^(0 + Pc(0 + PsiO + PpiO . (15) 
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In principle one has to diagonahze simultaneously momentum, energy and charge, i.e. 

P+\^i) = P+\^^) , P-I^-,) =P.-|^,) , and g3|^,)=0, (16) 

repectively. The denumerable eigenvalues of momentum and energy are denoted by and P~ , 
respectively. For the momentum and the charge this is not difficult, since any arbitrary Fock state 

W) = 4i42 • • • ^mi^ma • • • 4i42 ' • • 1°) (l*") 

with the same number of b- and d-particles is a solution, labelled appropriately by = 1, 2, . . . , A^. 
The number of states is finite for any fixed value of the 'harmonic resolution' Q il' = LP^ /tt. 
Since energy and momentum commute, they span a complete set of states for diagonalizing the 
energy. The Fock space vacuum |0), defined by 

an |0) = , 6„ |0) = , and |0) = , (18) 

has no particle content and thus vanishing momentum. 

Next to being an operator in Fock space the Hamiltonian is a Schrodinger operator with respect 
to the gauge field v or Generating a complete set of functions V'n(C) by solving a Schrodinger 
equation with a convenient single particle potential V{C), 

(-4ff'^ + v{0)M0 = U^nMO , (19) 
the solutions to Eq.([I^) must have a non-separable structure, i.e. 

oo N 

m = J2T.ci%MC)W), with c»^c«c». (20) 

n=0 u=l 

This statement holds in general, since the coefficient functions of the Fock-space operators in the 
Hamiltonian depend on in a complicated way. The physical vacuum \vac) that is the state with 
lowest energy of the full theory will thus acquire structure with respect to ^ as in Q, see also 
entirely due to a 'physical gauge' like 5_A+ = as opposed to the 'light-cone gauge' A+ = 0. But 
these structures will be disregarded in the rest of the paper where the vacuum is simple due to the 
model assumptions. 

4 The Model within the Model 

In the present work we approach the problem from the tail. We are interested in the particle sector, 
its spectral properties and wavefunctions. As a first step, in order to face a tractable problem, we 
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do not neglect but omit by hand ao, and in addition suppress the fluctuations in the gauge mode, 
i.e. 

ao = and = 5,,o <5(C) , (21) 

see also Ref. 0]. For C = the Hamiltonian becomes then 



= ^Fock ^ (0) + (0) + P^(0) . 



(22) 



One should emphasize that this model does not correspond to the light-cone gauge since we keep 
gvL = tt/2 as the static value around which the dynamic treatment would fluctuate. Last not least, 
we perform a Tamm-Dancoff truncation by restricting to the lowest, the 2-particle sectors. This 
'model within the model' will provide useful insight into the structure of the full solution. 

The 'contraction terms' Pq as given in Appendix ^ should deserve some explanatory remarks. 
In a first step one gets for them rather straightforwardly 



-P- 



E 

n=l 



^n{C) + Jmi+C) ,i , , Jmi — C) 



n 



d'^ d 



(23) 



The 'self-induced inertias' or 'tadpole diagrams' / and J contain logarithmically diverging pieces, 

2 4n 4n 



^(C) = /ig + 5'E 



+ 



+ 



n=l 
oo 



(m + C) ' (m-n + C)^ (m + n + C)^ 

2 4n 4n 

(m-C) {m-n-C,f {m + n-C,f 

2 ^ 4 (m -I- C) 4 (m + C) 

n [m — n + ^)^ (m + n + (^)^ 



(24) 



4(m + C) 4(m + C) 
(P-C) ' (P + C) ' {m-pf {m+pf 



+ 



+ 



(25) 



through terms like X^n^i ^/n. They can be cancelled by renormalizing the mass like 



^il = ^x''- IQg" - 2^2 



E - + E- 



1 



(16 + 67)52-4525]-, 



n=l 



m=l 



n 



(26) 



with 5 = g/^J\Q■K. Important is that the renormalized mass jjL is finite. Its value is not unique, and 
thus both the '16' and the constant 



00-1 00-1 
C^2 y --2V 



41n2 



m 



n=l 



n 



(27) 
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could be absorbed into the renormalized /x. But here is the problem: The sums in Eqs. (|2^ ) and 
(p5|) are not completely symmetric in the integers and half- integers. Converting the expressions 
tabulated in Table || generates the constant C in either one place or the other. In the continuum 
limit the two sums in Eq. ( |27| ) tend to cancel in line with Section |6[ We therefore shall work with 
C = whenever not mentioned otherwise. 

5 Numerical Solutions 

Consider two types of 2-particle Fock states, referred to as the 'aa'- and the 'bd'-space, respectively, 

K>a = 4a/f-nlO) foi^ n = l,..., 
and |m), = 6^4-^10) = ^> • • ■ > ^ " ^ • (28) 

They are orthogonal, {m\n) = 0, and simultaneous eigenstates of and P"*", i.e. 

Qz\n), = Q , P+|n), = ^i^|n)„ , (29) 
Q3|m), = , P+|m), = ^K|m), , (30) 

for all values of n and m, as they should according to Eq.(|T6|). Gauss' equation (^ prevents us 
from including also 'charged' states like for example ^m'^t-m I*-*)- Any of the K + [K/2\ linear 
superposition of these basis states like 

[fl ^ K-l ^ [f] K-l 

= E ("I l^)a + E ("^1 \i) \m)^, = E {n\^i,a) \n)^ + ^ iM^i.h) \m)^, (31) 

"=1 m=i ^1=1 m=i 

are thus eigenfunctions of and with the same eigenvalues, while the unitary matrix U = 
Ca + Ch determines itself by diagonalizing P~, with the {K + [K/2\) eigenvalues and subject 
to normalization 

{i\U^U\i) = l, foralH. (32) 

Since P'^Py = 2P^P^ is to be interpreted as the operator of invariant mass squared M^, the 
eigenvalues of P^ will be presented below as the product 2P^ P^ = Mf . As unit of mass we shall 
use mu = g = S'/VT^Gvr. Note that the coupling constant g has here the dimension of a mass. 

As a rather welcome advantage of a DLCQ calculation, the diagonalization of 100 x 100 matrices 
takes only a couple of milliseconds on a modern work station. One thus could go to fairly large 
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Figure 1: Invariant mass squared eigenvalues Mf versus Harmonic Resolution 
K. - Left side (a): The Fock space is restricted to only the aa-states \n)^. Note 
the almost complete degeneracy of all K/2 states. Right side (b): The Fock space 
is restricted to only the bd-states |n)^. - Note: All eigenvalues increase roughly 
linear in K, and show no trend to stabilize. - Parameters values are fj? = and 
ml = 2(100m„)2. 

values of the resolution within our simple model. However, in order to unravel the structure of the 
spectra and not to be overwhelmed by a flood of data, we first keep the mass parameter fixed to 
fi = and restrict ourselves to 20 < < 80. In view of typical values for the harmonic resolution 



K < 22 quoted in the literature [y^, |ll|, |l^, |lj] the above range is still large for practical purposes. 



1. The reduced 2-particle Fock spaces. When one restricts to only the aa-space, the Hamiltonian 
matrix (n'| P~ \n) is diagonal from the outset, since the a-particles have no interaction according to 
Table They only have self-induced inertias. To first order of approximation, for moderately small 
values of /x these inertias are linear in the momentum. The eigenvalues of P~ become approximately 
independent of K and thus the invariant masses therefore linear in K and highly degenerate, as 
displayed in Figure |l]a as function of K. 

If one restricts oneself to only the bd-space, the Hamiltonian matrix is non-diagonal due to the 
interaction between the b- and the d-particles. As displayed in Figure lb, this causes the eigenvalues 
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Figure 2: Invariant mass squared eigenvalues Mf versus Harmonic Resolution 
K. - The Fock space includes both the aa- and bd-states. Note: The lowest 
eigenvalues are roughly independent of K, with a roughly constant spacing. - 
Parameters values are /u^ = and = 2(100mu)^. 

to spread over a wider range of energies but obviously not strong enough so as to stabilize the low 
energy parts of the spectrum. One understands this behaviour analytically, see Section |6|. 

2. The 2-particle Fock space. The behaviour changes drastically if one includes simultaneously 
the aa- and the bd-states. One now allows for the virtual scattering from the aa-space into the 
bd-space and back, represented by the matrix elements Saa{ni,n2;ns,'n4) in Table ^. In Figure 
the first 19 mass squared eigenvalues are plotted as functions of K. At a first glance the eigenvalues 
are iC-independent, with a roughly equal spacing like the states in a harmonic oscillator well. At 
the left upper corner of the figure seem to develop some 'crossings' of the states, which we shall 
resolve and explain below. At a second glance, the spacing is not completely independent of the 
energy, rather it decreases slowly with increasing energy, typically like the eigenstates in a potential 
with linearly increasing walls. 
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We like to emphasize the discreteness of the spectrum. Discreteness surviving the continuum 
limit K ^ oo should be considered as the earmark of confinement, even if such a statement seems 
premature within the present 1+1 dimensional 'model within a model'. Most if not all field theories 



spectrum at all, in a model with no other ingredients than pure gauge theory. 

A word of caution seems in order. The restriction to the 2-particle sector is ad hoc and should 
not be overemphasized. Including more and more particles, one would expect structures which 
would tend to true continua for ever increasing K. At least parts of the so obtained spectra would 
correspond to 'physical glue-balls in relative motion'. In short and over-stressing the point, one 
would expect spectra which qualitatively resemble the multi-particle spectra of QED in (1+1) 
dimensions, as displayed in the figures of Ref.|jl5|, |l6|. 

3. Symmetries and Multiplet structure. In the naive thinking of color invariance a color singlet 
should have a finite spectrum. Are our eigenstates color singlets? Opposed to this interpretation is 
that the spectrum diverges when restricting to either aa- or bd-states alone as displayed in Figure 
Only when both are included, one has a chance for having wave functions which are invariant under 
rotations in color space. Can one be more specific? What are the symmetries of our problem? 

By inspecting the numerical results for the eigenfunctions one observes three distinctly different 
classes: (I) those with ^i^a ~ ^i,b, (H) those with ^i^a ~ 0, and (III) those with ^i^a ~ —'^^i,b, 
for comparable size of the single particle momenta. A separation into two groups can be based 
on an exact symmetry: the Hamiltonian is invariant when exchanging the b- and the d-particles. 
Denoting the corresponding charge-conjugation operator by Cjr, its eigenvalues must be = ±1. 
In fact all numerical eigenfunctions are either charge-conjugation even or odd, by inspection, and 
the latter coincides with the class II states. Aiming at a measure to classify the states we introduce 
the operators 



particularly QED confine in 1+1 dimension 1£]. Remarkable however is that one gets such a 



oo 



oo 




(33) 



Together with Q3 they satisfy SU(2) commutation relations in the Weyl representation 






(34) 



with the group invariant 



^Q3Q3 + Q-Q+ + Q+Q 



(35) 
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20 40 60 ,^ 80 20 40 60 80 
K K 

Figure 3: Spectral Decomposition with Multiplet Q. - For each state \i,K) 
the multiplet expectation value Q is calculated, and the spectrum disentangled 
correspondingly, see text. - Upper left (a): Full spectrum. Upper right (b): 
Spectrum of the quasi-singlets (Q = 0). Lower left (c): Spectrum of the quasi- 
triplets (Q = 1). Lower right (d): Spectrum of the quasi-pentuplets (Q = 2).- 
Note: Only the singlet states have a stable spectrum. - Parameters values are 
= and ml = 2(100m„)2. 

The eigenfunctions of the Hamiltonian cannot be eigenfunctions to since the two do not com- 
mute, see also Section]^ below. One can however calculate the expectation values qi = 
The 'effective eigenvalues' Q^, defined by QeiQe + 1) = Qi^ turn out to be close to the numbers 
0, 1, 2, . . .. We therefore define 

{0, if 0.0 < < 0.3 (class I); 
1, if 0.7 <Qe< 1.3 (class II); (36) 
2, if 1.7 < Qe < 2.3 (class III). 
It is remarkable, that none of the expectation values drops out of the comparably narrow limits 
set in these equations. Is this a consequence of a residual symmetry in the Hamiltonian? - For 
2 particles the largest possible value is Q = 2. As it turns out, all numerical eigenstates have 
a charge-conjugation parity = (—1)*^, particularly the triplet is charge-conjugation odd. In 
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Figure 4: Structure Functions versus Momentum Fraction x. - The structure 
function is calculated both for the a- and the b-particles, (x = n/K, squares) 
and (x = m/K, bullets), respectively, for the ground and the first excited state. - 
Note: The curve interpolating between squares and bullets is remarkably smooth. 
Note also, that the structure function of the excited state has more nodes than 
the one of the ground state, in accord with expectation. - Parameters values are 
Ij"^ = andK = 60. 

Figure ^ the same spectrum as in Figure ^ is displayed, but separated according to the multiplet- 
value Q. As it turns out, the crossing situation mentioned in the context of Figure Q is due to the 
crossing between Q = and Q = 1 states. All eigenvalues of the singlet Q = become virtually 
independent of the resolution. As opposed to this all members of the triplets or pentuplets, Q = 1 
or Q = 2, respectively, are at least linear in K, tending to infinity in the continuum limit. 

The latter result is remarkable since it is one of the few direct evidences that only singlets can 
have finite masses in a non-abelian theory. Some subtleties related to this interpretation will be 
discussed further below in Section |^. 

4. Structure Functions. Next to the spectrum we investigate explicitly the eigenfunctions. Due 
to diagonalization one knows them for all eigenstates, particularly their projections on the Hilbert 
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space, {n\^i) = (n| Ca \i) and {m\^i) = (m| Cb see also Eq.(^). The probabilities to find an a- 
or a b-particle with longitudinal momentum = "-j = rrij- are given by 

alan = {n\ Ca \if , and blhm \^^) = {m\ Cb \if , (37) 

repectively, and are related to the structure functions f{x) by 

ft\n/K) = K{^i\ alan , and ft\rn/K) = K blb„, . (38) 

In the continuum, the probability to find a parton with longitudinal momentum fraction between 
X = n/k and x + dx is dx f{x). Note that the structure functions are not normalized to unity, 
because of Eq.(|32[), rather they obey the sum rule Jq dx (^^fj^"'\x) + fi''\x)j = 1. In displaying 
them in Figure ^, we restrict ourselves to the first two states. Due to the residual SU(2) symmetry 
the structure functions for the a- and b-particles turn out to be extremely similar, for which reason 
they have been compiled in the same plot. By inspection the structure function for the first state 
can be approximated by 



u^) = o .± "2 [^(i-^)r , (39) 



2 r(2 + 2a ) 
3[r(l + a) 

with an exponent numerically closer to a ~ | than to a ~ |. All higher states have more nodal 
structures. 

5. Mass-Dependence. The results presented thus far have been calculated for a vanishing value 
of the mass parameter ^u. For a mass parameter /i much larger than the effective coupling constant 
g one can omit the interaction to first order of approximation. With increasing /i, the spectrum 

2 

must become more and more like the one of two free massive bosons, that is like ~ —tt — v. Cor- 

' x(l—x) 

respondingly, the structure function of the lowest state will be peaked at x = 1/2 with an ever 
decreasing width like a 5-function. As shown in Figure ^, this is exactly what happens. Keeping 
g fixed and increasing /i the spectrum changes from being 'interaction dominated' at /i ~ in the 
manner as having been shown in Figure ^ to being 'mass dominated' beyond fi ~ WOniu- A good 
earmark of the latter is the decreasing level density with increasing excitation. 

6. Does the Continuum Limit exist? The spectra shown in Figs. ^ and ^ for 20 < X < 80 did 
appear to be stable as function of K. Are they really? What happens when one increases K 
by order of magnitudes? Does the continuum limit exist? - The question can be answered either 
analytically by converting the Hamiltonian matrix equation to an integral equation, or numerically 
by increasing the order of magnitude of K. We have done both. The continuum limit is discussed 
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Figure 5: Invariant mass squared eigenvalues Ei = 2Mj^/(1000m^ + 40^^) versus 
Mass Parameter fj,'^. Note the logarithmic scale and how the spectrum changes 
from being 'interaction dominated' at smaller to being 'mass dominated' at larger 
values of ii. - The harmonic resolution is fixed at K = 50. 

in detail in Section Here we present the numerical results, by varying up to ~ 1000. 
Plotting the results in Figure |6| on a logarithmic scale one is able to unravel logarithmically small 
variations which skip the observation on a linear scale. The unexpectedly large but still modest 
variation of the eigenvalues with increasing resolution particularly their positive curvature seems 
to indicate convergence to the continuum limit, in line with the analytic considerations below in 
Section ^ We have verified this behaviour numerically also for larger values of the mass parameter 
^, but renounce to display the results. For negative values of the mass parameter, however, i.e. 
for /i^ < the eigenvalue curves tend to develop an increasingly strong negative curvature with 
respect to InK, again in line with analytic considerations. 
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Figure 6: Eigenvalues Ei versus Large Harmonic Resolution K. - The lowest 
three eigenvalues are plotted in however different mass units Ei = P~ /mf 
with ml = 75m^, m2 = 480m^, = 730m^. They are reasonably but not 
completely stable and vary on a level of less than 10 per cent. They seem to 
converge in the continuum limit X — > oo. 

6 Coupled Integral Equations in the Continuum Limit 



According to Eg. (pq) , one deals with orthogonal vector spaces {m\n) = 0. In such a space the 
eigenvalue problem P^ \ = becomes a set of two coupled matrix equations, i.e. 



K- 



^ {n\p-\m'){m'\^i) = [pr - {n\ p- \n)^ {n\^ i) , (40) 



[K/2] 



^(m|P (n'l^'i) + ^ (m| P~ |m') (m'|^'j) = (p." - (m| |m)) (m|^'i) , (41) 



n'=l 



2 



see also Eq.(pll). As mentioned, matrix elements like (n| P~ |n') vanish in the present model. Hence 
forward we substitute P^ = uJi/P^^ and drop the index i. The notation accounts for the obvious 
fact that the off-diagonal matrix elements are calculated differently from the diagonal ones and 
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that the latter are taken to the r.h.s. of the equations. Correspondingly, in the continuum limit, 
one is confronted with a set of two coupled integral equations 



UJ 



^a{x) = Ia{x) + -f^dx' {x\ Ka \x') V'b(x') , (42) 

JO 

= h{x) Mx) + i^dx' {x\ Kb \x') Mx') + -f'dx' {x\ Ka \x') Mx') ■ (43) 
Jo Jo 



More specifically, the continuum limit is obtained by the limiting procedure n ^ oo and K ^ oo 
at fixed momentum fractions x = n/K, replacing sums by integrals like for example 



(44) 



2 .1 

F{m,m') = -f dx' F{x,x') , with F{x,x') = K F{m,m') . 

Care has to be taken of a proper removal of the 'diagonal kernels', i.e. 

-l^dx' F{x, x') = r 'dx' F{x, x') + dx' F{x, x') , (45) 

Jo Jo Jx+€ 

with the limit e — > taken at the end. Matrix elements are thus related to kernels by 

{x\Ka\x') = KP+y {n\p- \m') , and {x\ Kb\x') = K P'^y {m\ p- \m') . (46) 

From Table || one gets 

{x\K W) = {^x + x'){2-x-x') 2g^ + x - x') {I - x + x') 

"'^^ [x-x'f ^x{l-x)x'{l-x') [l-x-x'f ^x{l-x)x'(l-x') ' 

{x\K\x') = {x + x'){2-x- x') (1 - 2x) (1 - 2x0 

^'^^ {x-x'f ^x{l-x)x'{l-x') ^x{l-x)x'{l-x') ■ ^ ' 

In the sequel the obvious symmetry {x\ Ka \x') = {x\ Ka |1 — x') will be used without further men- 
tioning. Note the non-integrable quadratic singularities at x' = x and x' = 1— x. The wavefunctions 
become {m\^) = ^a(x) and {n\^) = V'b(x). The diagonal terms are obtained from Table [l|, i.e. 

We now seek the solution in the regime which corresponds to the color singlets, namely 



tpb{x) = -06(1 - X) = Tpaix) = 1p{x) . (50) 

The coupled integral equations (|^) and (|3|) degenerate then into two equations for one function 



u)'4>{x) = Ia{x) ip{x) + + dx {x\ Ka\x') Tpix') , and (51) 

JO 

u) V'(x) = Ib{x) ipix) + / dx' (x| Kb \x') V'(x') + --[ dx' (x| Ka \x') V^(x') . (52) 

Jo 2 Jo 
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One of them must therefore be redundant. To solve a singular equation like (|5^) one adds and 
subtracts a term like I(x) = Jq dx' {x\ Ka \x'), i.e. 



LO l\j{x) = {I{x) + Ia{x)) tp^x) + -f dx' {x\ Ka \x') {lp{x') — 1p{x)) 

Jo 



(53) 



Close to the singularity, ip{x) — ip{x') ~ x — x' converts the quadratic singularity into a pole which 
can be integrated by some principal value prescription. Now look at 



T(\^T(\ [l + l6^±C^^,.2f'dx' 
I{x) + Ia{x) = — + Ag 7;— ;:,^2 



(x + x') {2 — x — x') 



\/x{l — x)yjx'{l — X') 



(54) 



x(l — x) ' ^ Jo — x'Y 
Setting in the square bracket x' = x + A and expanding with A one notes that the first non- 
vanishing term is oc A^. This cancels the singularity, and the integral becomes finite. The integral 



equation (51) has thus a solution. - As for Eq.(|52D we subtract the latter from the former to get 



Cg' 



X 1 



■ ^(x) = -f dx' ip{x') 
Jo 



{x\Kb\x')--{x\Ka\x') 



(55) 



Inserting the kernels from Eqs. ( |47| ) and ( |4q ) the singular terms cancel precisely. One remains with 

Cg^ , , o (l-2x) . , (l-2x') , , 

2.= ^^//x>(x')-i_i==0. (56) 



■ ip{x) 



x(l-x)'''"' a/x(1 - x) Jo \/x^) 

because the integrand is an odd function around x' = ^. One has to conclude: The second integral 
equation is consistent with the first one iff C = 0. Otherwise the continuum limit does not exist. 
This is in line also with the statement from color invariance: All three particles should have the 
same mass, see Eg. (^9[) . - One remains with Eq.(^) which units adjusted is identical with the one 
of Bardeen et al. jl^ and of Klebanov et al. [|lT|, although the latter had been derived in another, 
i.e. in the light-cone gauge. We refer to their work particularly with respect to the endpoint 
analysis of the solutions. 



7 Charge Conservation and Color Singlets 

We take the multiplet structure of the spectra in Section |5| as kind of an empirical fact, and wonder 
whether they can be put on more solid grounds. This turns out more difficult than anticipated. - 
Of course, one always can associate charges with the currents as defined in Section |5|, particularly 

Q- = J dx^Jl, = Q^_, and the familiar Q3. Upon evaluation they become 

00 00 
Qs = ^ dl^dm - bl^bm , and Q- = AoUobi) + ^ Analb^_^_i - Bnand^^_^ , (57) 
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with {aobi)s = ^{aobi + biao). In the fundamental modular domain the coefficients are 



n 



, and Bn = ]^i \ - — - + 

2 \ \ n \ n — z 



(58) 



Since 



O'O , Qz = 0, see Ref.|§], one calculates the SU(2)-commutators as follows: 



Qs , Q- 



Q- and 



Q- , Qh 



Q3 = A(z,ao) 



(59) 



The right hand side of this equation turns out to be 
4 



-A(z,ao) = 



aobi) , (aob\) 

2/s \ ^Js 



+ 



+ 



+ E 



- z2 V "~ 2 



n=l 



n=l 



+ + Onfln ) (60) 



2 



-, ~ X z - z -/ n (n^ — z2) \^ 

n=l ^ ' n=l ^ ' ^ 

In passing one notes that the charges of Eq.(|3^) agree with those of Eq.(^7|) for z = (thus 
A = 0), and for oq = 0. Can one anticipate for the general case that A{z, oq) = 0? - One has two 
perspectives for the future: Either (I) the explicit solution for ag renders A = 0, or (II) one requires 
A = as a subsidiary condition for solving oq- Here the matter rests as long as one cannot solve 
explicitly for oq. 

But in either case seems to be a problem: The charges as defined by Eq.(|57|) cannot be conserved 
since df^J^ ^ 0. True, our model has a conserved four-current df^J^ = 0, since this follows directly 
from the color-Maxwell equations (9qF"^ = gJ^. But the two currents are not identical since 
J/5 = J/3 + 1 [F/5° , Aa]. Only for the latter holds 



d 



+L 



Q = 0, with Q 



(61) 



Thus Q = Q + 2iL[d+A+ , A+] + 2Lg[[A+ , (A-)^] , A+], or 

Q3 = Q3, and Q_ =Q_ -2Lc/i;2(^;)g, 



(62) 



But now, in the process of inverting the Gauss' equations (|9|), their zero modes have been required 
to vanish [^. This in turn leads to Q± = 0, in total and unpleasant opposition to tribal beliefs 
that the true and dynamically conserved charges Qa obey a non-Abelian group structure. 
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As a possible way out we proffer a different gauge choice namely 

(A-)o = 0, (63) 

in addition to Eq.(^). We leave it for future work to check whether this choice is possibly in conflict 
with the first gauge condition in Eq.(|3|). But the proffered choice would have the advantage to 
be manifestly invariant under color rotations and to treat the zero mode components of A~ the 
same way in all three Gauss' equations (|9|). Moreover, since Q = Q, both charges would be strict 
constants of the motion, reconciling thus naive and refined considerations. 

But perhaps these considerations are academic since all zero modes become sets of measure zero 
in the continuum limit. The solutions of the integral equation, for example, are color singlets in 
the strict sense. As mentioned, they are independent of the gauge choices in the zero mode sector. 



8 Summary and Discussion 

The approach taken here follows closely some earlier work Q. Beginning with SU(2) pure gauge 
theory in (2+1) dimensions in the front form one suppresses the transverse coordinate dependence 
of the gluons and obtains a (1+1) dimensional gauge theory coupled to adjoint scalar matter. The 
present work has one major aim: We would like to get a first and rough idea on how the excitation 
spectrum as well as how the structure functions look in such a dimensionally reduced quantum field 
theory. Therefore, in order to have a tractable formalism and in contrast to earlier Q and ongoing 
work we suppress here by hand the dynamical impact of the topological gauge zero mode ((") 
and of the constrained zero mode (ao)- We perform furthermore a Tamm-Dancoff truncation and 
include only the 2-particle Fock space. 

What are the results? Most important we think is the result displayed in Figure ^, namely that 
the solutions corresponding to color singlets {Q = 0) have a discrete and finite spectrum over a 
wide range of the harmonic resolution K. The states corresponding to color triplets (Q = 1) or 
pentuplets {Q = 2) can well be separated and tend to have a very large if not infinite mass. This 
appears to be one of the few available concrete pieces of evidences that only the color singlets can 
have finite mass in a non-abelian field theory. The method to identify these states is somewhat 
pragmatic since the model assumption on the zero mode, oq = 0, violates strict SU(2) invariance. 

In Section ^ the continuum limit was established by deriving an integral equation in the manner 
of Bergknoff [18| for QEDi+i, see also |15]. As demonstrated, the continuum limit exists only for 
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C = 0, resulting for the color singlets in the integral equation 

,, , + 16g^ ,, , ,^of^ dx' . , , ix + x') [2 — X — x') , , , 

uJi^{x) = ^ ^^{x)+4gU- ^ 4V X - \ Mx') . 64 

x{l — x) Jo{x — x'y — x) x'(l — x'j 

The integral equation is identical with the one derived earlier by Bardeen et al. |0] and by Klebanov 



et al. 1 11], a remarkable fact in view that both these authors have used the light-cone gauge as 
opposed to the present light-cone Coulomb gauge. 

The present model is restricted to the lowest non-trivial Fock states, those with particle number 
two. Fock states with three partons like ajj6jjjd|. |0) are ruled out since all matrix elements between 
the 2- and the 3-particle sector vanish due to color invariance. The next higher Fock-states have 
thus four partons: aaaa-, aabd, and bdbd-states. Based on the numerical experience with the 



massive and massless Schwinger model |15, 16 1, one can guess that the admixtures of the 4-parton 
states in the low lying part of the spectrum are of the order of 10^'^ and thus probably negligible. 

A rather interesting feature of the model should be mentioned. For sufficiently small variations 
of the mass parameter fi around the value zero the spectrum as displayed in Figure ^ can be moved 
up and down unchanged almost at gusto, including the case that the lowest eigenvalue coincides 
with zero. A trace of that survives even in the larger variations of Figure ^. It looks as if a mass 
difference like the one in Eq. ( p^ ) is a dimensionless number characteristic for the present model. It 
should be interesting to get an analytic estimate, because the mechanism allows to generate huge 
mass ratios M2/M1 as they are characteristic for hadronic physics. 

What can one learn from this work beyond our particular play model? Let us return to the 
fundamental assumption 9jA^ = as in Q. Dimensionful quantities are never strictly zero, so let 
us be more quantitative. In the full 3-|-l dimensional treatment one would introduce an invariant 
mass cut-off like for example the one of Lepage and Brodsky, see Q , 

It is covariant, the sum runs over all partons v. The scale A is at our disposal. For two massless 
gluons one obtains thus k\ < x(l — x)A^. Maximizing this by x = 1/2 and using transverse periodic 
boundary conditions with k = nvr/L^ one gets |n| < AL±/{2tt). Therefore, all transverse momenta 
are cut out when one defines the transverse length by 

L±^l- (66) 
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Only the transverse zero modes survive, and this is precisely the present model. 

The fact that we have started from 2+1 dimensions weighs less. In 3+1 dimensions, the model 
with only transversal zero modes in the manner of 0] is only marginally more complicated due to 
the non-abelian commutator term. First estimates show that its impact is not dominant. With 
a grain of salt, the present color singlet solutions would then correspond to the helicity aligned 
2++-glue-balls. If one sets the cut-off scale as the typical hadronic energy A ~ 1 GeV one gets for 
the transversal size L± ~ 0.6 /m, a value not untypical for hadronic sizes. 

One can push these simple considerations even further. The present coupling constant g is 
related to in 3+1 dimensions |0] by gs = 2gL±. Our unit of mass niu = g/A^^Tr can thus be 
expressed in terms of = g'^/^ir and A. A value of ~ 0.4 looks reasonable as compared to the 
empirical as(Mz) — 0.12. With Eq.(|6^) one obtains 

rriu = V^s^ ~ 25 MeV , (67) 

which allows to convert our numerical results to physical units. By order of magnitude, one reads 
off Figure | for the first two states ~ 0.0018 and M| ~ 0.038 mg, thus 

-\ fM| - M?) ~ 800 , Ml ~ 150 MeV , M2 ~ 600 MeV . (68) 

These numbers are not untypical for hadrons. 

In conclusion, the physical picture emerging from this work might over-stress the point but 
without being necessarily false: Hadrons particularly glue-balls have no transversal structure up 
to transversal sizes of about 0.5 fm, compare also with the recent work of van Baal ||T^. The 
field lines are parallel due to the basic ingredient of the model, in accord with phenomenological 
flux-tube models. Up to these sizes, the structure resides in the longitudinal direction as mirrored 
in the structure functions displayed in Figure |^. 
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A The Matrix Elements of the Contraction Energy 



Tabic 1: The contraction part Pq of the Hamiltonian is expressed in terms 
of Fock-space operators. The couphng constant is absorbed into the coefEcient 
= — — . As discussed in the text, one should use C = in the calculations. 

lOTT 



Graph 


Matrix Element 


n \/\/\/\/\/V. 




^ m 


.UC) = ^P + ^2(^i6 + 8(TO + ofE (^^ + '"e ^)+4(f^ 

+ C E „(„2_^2) 1 C(™ 1 C) ^ (k^-c^Y ) 

Ti= ^ A; — m + 1 / 
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B The Matrix Elements of the Seagull Energy 



Table 2: The seagull part Pg of the Hamiltonian is expressed in terms of Fock- 
space operators. The coupling constant is absorbed into the coefficient = —— . 

lOTT 

The indices of a and are integers, those of b,b\d and S half-integers. The 
Kronecker delta (^^g+JJ^ referring to momentum conservation in S{ni,n2;n3,n4) 
is not kept track of exphcitly. 



Graph 



Matrix Element 



"l+C I / "3+C 

"3+c V "1+c 



n2-C 



+ 



rii-C y n2-( J (rai-ri3)2 




+ 



"2-C 

ni+C 



ni+C 

"2-C 



n3+C 



1 kAAAA/V 



(rii-ri3-C)^ 



+ 



"2 + C 



»2+C 

"1 



n4 

m+C 



2g' 



ns+C 



IvA/WV 



/ rij-C ] i-2g^) 



-P- 



C) ^r!2 ^"3 

bns dj^^ + d^^^ 6tyi2 ^ni ) 
+ E Sabini,n2;n3,n4;C) al^bl^bn^a^^ + Sabini,n2;n3,n4; -() al^dl^dn^a 
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C The Matrix Elements of the Fork Energy 



Table 3: The fork part Pp of the Hamiltonian is expressed in terms of Fock-space 



operators. The coupling constant is absorbed into the coefficient g 
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■. The 



indices of a and are integers, those of h,h\d and d'^ half-integers. The Kro- 
necker delta (5JJ2+n3+„4 referrf 
is not kept track of exphcitly. 



necker delta 6'^^_^^^^^^ referring to momentum conservation in -F(ni; n2, ns, 7x4) 



Graph 



Matrix Element 




Fba{ni; n2, ris, n4] C) = 



ni+C 
"2+C 



+ 



»2+C 



rn-C 



2£_ 



n3 + C J (ni-n2Y 



,4 



»i+C 



+ 



"2 



n^+C _ 

ns \l 114 +C ) (ni-n2-\-CY 



Fab(ni; n2, n3, C) 



+ 



»2+C 



»3-C 



+ 



+ 



n2 + C 



n3-C 



n4 

"3-C 



"2+C 



"3-C 

114 / (ill -112 -C)^ 



»2+C 

"4 / (ni-n^+C)^ 



2g' 



F^^ 



-Pit = J2Fba(ni; 112,113, n4;() bl^ 

Fit 

^713 '^rii + Fbbiiii; 112, 113, 114, -C) c^ij 

Fit 

+ J2Fabini;n2,n3,n4; () al^bn^dnsam + h.c. 
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